Footprinting molecular electrostatic potential surfaces for calculation of solvation energies.

A liquid is composed of an ensemble of molecules that populate a large number of different states, so calculation of the solvation energy of a molecule in solution requires a method for summing the interactions with the environment over all of these states. The surface site interaction model for the properties of liquids at equilibrium (SSIMPLE) simplifies the surface of a molecule to a discrete number of specific interaction sites (SSIPs). The thermodynamic properties of these interaction sites can be characterised experimentally, for example, through measurement of association constants for the formation of simple complexes that feature a single H-bonding interaction. Correlation of experimentally determined solution phase H-bond parameters with gas phase ab initio calculations of maxima and minima on molecular electrostatic potential surfaces (MEPS) provides a method for converting gas phase calculations on isolated molecules to parameters that can be used to estimate solution phase interaction free energies. This approach has been generalised using a footprinting technique that converts an MEPS into a discrete set of SSIPs (each described by a polar interaction parameter, εi). These SSIPs represent the molecular recognition properties of the entire surface of the molecule. For example, water is described by four SSIPs, two H-bond donor sites and two H-bond acceptor sites. A liquid mixture is described as an ensemble of SSIPs that represent the components of the mixture at appropriate concentrations. Individual SSIPs are assumed to be independent, so speciation of SSIP contacts can be calculated based on properties of the individual SSIP interactions, which are given by the sum of a polar (εiεj) and a non-polar (E(vdW)) interaction term. Results are presented for calculation the free energies of transfer of a range of organic molecules from the pure liquid into water, from the pure liquid into n-hexadecane, from n-hexadecane into water, from n-octanol into water, and for the transfer of water from pure water into a range of organic liquids. The agreement with experiment is accurate to within 1.6-3.9 kJ mol(-1) root mean square difference, which suggests that the SSIMPLE approach is a promising method for estimation of solvation energies in more complex systems.

Interplay of self-association and solvation in polar liquids.

The association constants for formation a 1:1 complex between 4-phenyl azophenol and tri-n-butylphosphine oxide were measured in mixtures of n-octane and n-decanol, n-octane and n-hexanoic acid, and n-octane and 2-ethylhexyl acetamide. The experiments provide insight into the competition between solvent self-association and solvent-solute interactions in these systems. The solvation properties of the three polar solvents are quite different from one another and from polar solvents that do not self-associate. Carboxylic acids form dimers in concentrated solution (>1 mM in alkanes). Carboxylic acid dimers have exposed H-bond acceptor sites that solvate H-bond donor solutes with a similar binding affinity to carboxylic acid monomers. The carboxylic acid H-bond donor site is inaccessible in the dimer and is not available to solvate H-bond acceptor solutes. The result is that solvation of H-bond acceptor solutes is in competition with solvent dimerization, whereas solvation of H-bond donor solutes is not. Secondary amides form linear polymers in concentrated solution (>10 mM in alkanes). The solvation properties of the secondary amide aggregates are similar to those of carboxylic acid dimers. Solvation of H-bond acceptor solutes must compete with solvent self-association, because the amide H-bond donor site is not accessible in the middle of a polymeric aggregate. However, the amide aggregates have exposed H-bond acceptor sites, which solvate H-bond donor solutes with similar binding affinity to amide monomers. Alcohols form cyclic tetramers at concentrations of 100 mM in alkanes, and these cyclic aggregates are in equilibrium with linear polymeric aggregates at concentrations above 1 M. The alcohol aggregates have exposed H-bond acceptor sites that solvate H-bond donor solutes with similar binding affinity to alcohol monomers. Although the alcohol H-bond donor sites are involved in H-bond interactions with other alcohols in the aggregates, these sites are sufficiently exposed to form a second bifurcated H-bond with H-bond acceptor solutes, and these interactions have a similar binding affinity to alcohol monomers. The result is that self-association of alcohols does not compete with solvation of solutes, and alcohols are significantly more polar solvents than expected based on the properties of alcohol monomers.

Influence of solvent polarity on preferential solvation of molecular recognition probes in solvent mixtures.

The association constants for formation of 1:1 complexes between a H-bond acceptor, tri-n-butylphosphine oxide, and a H-bond donor, 4-phenylazophenol, have been measured in a range of different solvent mixtures. Binary mixtures of n-octane and a more polar solvent (ether, ester, ketone, nitrile, sulfoxide, tertiary amide, and halogenated and aromatic solvents) have been investigated. Similar behavior was observed in all cases. When the concentration of the more polar solvent is low, the association constant is identical to that observed in pure n-octane. Once a threshold concentration of the more polar solvent in reached, the logarithm of the association constant decreases in direct proportion to the logarithm of the concentration of the more polar solvent. This indicates that one of the two solutes is preferentially solvated by the more polar solvent, and it is competition with this solvation equilibrium that determines the observed association constant. The concentration of the more polar solvent at which the onset of preferential solvation takes place depends on solvent polarity: 700 mM for toluene, 60 mM for 1,1,2,2-tetrachloroethane, 20 mM for the ether, ester, ketone, and nitrile, 0.2 mM for the tertiary amide, and 0.1 mM for the sulfoxide solvents. The results can be explained by a simple model that considers only pairwise interactions between specific sites on the surfaces of the solutes and solvents, which implies that the bulk properties of the solvent have little impact on solvation thermodynamics.

Molecular recognition probes of solvation thermodynamics in solvent mixtures.

High-throughput UV-Vis experiments using four molecular recognition-based probes, made by the combination of two hydrogen bond acceptors, tri-n-butylphosphine oxide and N,N'-bis(2-ethylhexyl)acetamide, and two hydrogen bond donors, 4-phenylazophenol and 4-nitrophenol, were performed. The association constants for the 1 : 1 H-bond interaction involved in each probe system were measured in mixtures of a polar and non-polar solvent, di-n-hexyl ether and n-octane, respectively. Similar behaviour was observed for all four systems. When the concentration of the polar solvent was low, the association constant was identical to that observed in pure n-octane. However, once the concentration of the polar solvent exceeded a threshold, the association constant decreased linearly with the concentration of di-n-hexyl ether. Selective solvation in mixtures can be understood based on the competition between the multiple competing equilibria in the system. In this case, solvation thermodynamics are dominated by competition of the ether for solvation of H-bond donors. For the more polar solute, 4-nitrophenol, the selective solvation starts at lower concentrations of the polar solvent compared with the less polar solute, 4-phenylazophenol. Thus the speciation and hence the properties of systems containing multiple solutes and multiple solvents can be estimated from the H-bond properties and the concentrations of the individual functional groups.

The design, synthesis and pharmacological characterization of novel ß2-adrenoceptor antagonists.

BACKGROUND AND PURPOSE: Selective and potent antagonists for the ß(2) -adrenoceptor are potentially interesting as experimental and clinical tools, and we sought to identify novel ligands with this pharmacology. EXPERIMENTAL APPROACH: A range of pharmacological assays was used to assess potency, affinity, selectivity (ß(2) -adrenoceptor vs. ß(1) -adrenoceptor) and efficacy. KEY RESULTS: Ten novel compounds were identified but none had as high affinity as the prototypical ß(2) -adrenoceptor blocker ICI-118,551, although one of the novel compounds was more selective for ß(2) -adrenoceptors. Most of the ligands were inverse agonists for ß(2) -adrenoceptor-cAMP signalling, although one (5217377) was a partial agonist and another a neutral antagonist (7929193). None of the ligands were efficacious with regard to ß(2) -adrenoceptor-ß-arrestin signalling. The (2S,3S) enantiomers were identified as the most active, although unusually the racemates were the most selective for the ß(2) -adrenoceptors. This was taken as evidence for some unusual enantiospecific behaviour. CONCLUSIONS AND IMPLICATIONS: In terms of improving on the pharmacology of the ligand ICI-118,551, one of the compounds was more selective (racemic JB-175), while one was a neutral antagonist (7929193), although none had as high an affinity. The results substantiate the notion that ß-blockers do more than simply inhibit receptor activation, and differences between the ligands could provide useful tools to investigate receptor biology.

FieldScreen: virtual screening using molecular fields. Application to the DUD data set.

FieldScreen, a ligand-based Virtual Screening (VS) method, is described. Its use of 3D molecular fields makes it particularly suitable for scaffold hopping, and we have rigorously validated it for this purpose using a clustered version of the Directory of Useful Decoys (DUD). Using thirteen pharmaceutically relevant targets, we demonstrate that FieldScreen produces superior early chemotype enrichments, compared to DOCK. Additionally, hits retrieved by FieldScreen are consistently lower in molecular weight than those retrieved by docking. Where no X-ray protein structures are available, FieldScreen searches are more robust than docking into homology models or apo structures.

Thrombogenic collagen-mimetic peptides: Self-assembly of triple helix-based fibrils driven by hydrophobic interactions.

Collagens are integral structural proteins in animal tissues and play key functional roles in cellular modulation. We sought to discover collagen model peptides (CMPs) that would form triple helices and self-assemble into supramolecular fibrils exhibiting collagen-like biological activity without preorganizing the peptide chains by covalent linkages. This challenging objective was accomplished by placing aromatic groups on the ends of a representative 30-mer CMP, (GPO)(10), as with l-phenylalanine and l-pentafluorophenylalanine in 32-mer 1a. Computational studies on homologous 29-mers 1a'-d' (one less GPO), as pairs of triple helices interacting head-to-tail, yielded stabilization energies in the order 1a' > 1b' > 1c' > 1d', supporting the hypothesis that hydrophobic aromatic groups can drive CMP self-assembly. Peptides 1a-d were studied comparatively relative to structural properties and ability to stimulate human platelets. Although each 32-mer formed stable triple helices (CD) spectroscopy, only 1a and 1b self-assembled into micrometer-scale fibrils. Light microscopy images for 1a depicted long collagen-like fibrils, whereas images for 1d did not. Atomic force microscopy topographical images indicated that 1a and 1b self-organize into microfibrillar species, whereas 1c and 1d do not. Peptides 1a and 1b induced the aggregation of human blood platelets with a potency similar to type I collagen, whereas 1c was much less effective, and 1d was inactive (EC(50) potency: 1a/1b >> 1c > 1d). Thus, 1a and 1b spontaneously self-assemble into thrombogenic collagen-mimetic materials because of hydrophobic aromatic interactions provided by the special end-groups. These findings have important implications for the design of biofunctional CMPs.

Substituent effects on aromatic stacking interactions.

Synthetic supramolecular zipper complexes have been used to quantify substituent effects on the free energies of aromatic stacking interactions. The conformational properties of the complexes have been characterised using NMR spectroscopy in CDCl(3), and by comparison with the solid state structures of model compounds. The structural similarity of the complexes makes it possible to apply the double mutant cycle method to evaluate the magnitudes of 24 different aromatic stacking interactions. The major trends in the interaction energy can be rationalised using a simple model based on electrostatic interactions between the pi-faces of the two aromatic rings. However, electrostatic interactions between the substituents of one ring and the pi-face of the other make an additional contribution, due to the slight offset in the stacking geometry. This property makes aromatic stacking interactions particularly sensitive to changes in orientation as well as the nature and location of substituents.

Quantification of functional group interactions in transition states.

A new version of the double mutant cycle approach has been used for the evaluation of weak noncovalent interactions in transition states.

The role of the counteranion in the cation-pi interaction.

Chemical double mutant cycles have been used to quantify cation-pi interactions in chloroform as a function of the nature of the counteranion. The cation-pi interaction is -2.5 +/- 0.4 kJ mol(-1) and independent of the anion, even though the overall stability of the complexes varies by an order of magnitude due to competition of the anion for alternative binding sites.

An evaluation of force-field treatments of aromatic interactions.

Experimental measurements of edge-to-face aromatic interactions have been used to test a series of molecular mechanics force fields. The experimental data were determined for a range of differently substituted aromatic rings using chemical double mutant cycles on hydrogen-bonded zipper complexes. These complexes were truncated for the purposes of the molecular mechanics calculations so that problems of conformational searching and the optimisation of large structures could be avoided. Double-mutant cycles were then carried out in silico using these truncated systems. Comparison of the experimental aromatic interaction energies and the X-ray crystal structures of these truncated complexes with the calculated data show that conventional molecular mechanics force fields (MM2, MM3, AMBER and OPLS) do not perform well. However, the XED force field which explicitly represents electron anisotropy as an expansion of point charges around each atom reproduces the trends in interaction energy and the three-dimensional structures exceedingly well. Collapsing the XED charges onto atom centres or the use of semi-empirical atom-centred charges within the XED force field gives poor results. Thus the success of XED is not related to the methods used to assign the atomic charge distribution but can be directly attributed to the use of off-atom centre charges.

Substituent effects on cation-pi interactions: a quantitative study.

A synthetic supramolecular complex has been adapted to quantify cation-pi interactions in chloroform by using chemical double-mutant cycles. The interaction of a pyridinium cation with the pi-face of an aromatic ring is found to be very sensitive to the pi-electron density. Electron-donating substituents lead to a strong attractive interaction (-8 kJ/mol(-1)), but electron-withdrawing groups lead to a repulsive interaction (+2 kJ/mol(-1)).

Noncovalent Assembly of [2]Rotaxane Architectures.

Reversible zinc-pyridine coordination and hydrogen-bonding interactions have been used to assemble a [2]rotaxane from three components. Cooperativity in the macrocyclization process that results in the porphyrin dimer makes the system exceptionally stable. However, the kinetic lability of the zinc-porphyrin interaction means the dimer is in dynamic equilibrium with its monomer, and this has been exploited in the construction of a [2]rotaxane.